Kondo behavior, ferromagnetic correlations, and crystal fields in the heavy Fermion 

compounds CesX ( X=In, Sn) 
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We report measurements of inelastic neutron scattering, magnetic susceptibility, magnetization, 
and the magnetic field dependence of the specific heat for the heavy Fermion compounds Ceain and 
CesSn. The neutron scattering results show that the excited crystal field levels have energies Ei = 
13.2 meV, E2 = 44.8 meV for Cegln and Ei = 18.5 meV, E2 = 36.1 meV for CegSn. The Kondo 
temperature deduced from the quasielastic linewidth is 17 K for Cesin and 40 K for CeaSn. The 
low temperature behavior of the specific heat, magnetization, and susceptibility can not be well- 
described by J=l/2 Kondo physics alone, but require calculations that include contributions from 
the Kondo effect, broadened crystal fields, and ferromagnetic correlations, all of which are known 
to be important in these compounds. We find that in Ceain the ferromagnetic fiuctuation makes 
a 10-15 % contribution to the ground state doublet entropy and magnetization. The large specific 
heat coefficient 7 in this heavy fermion system thus arises more from the ferromagnetic correlations 
than from the Kondo behavior. 



I. INTRODUCTION 

In heavy fermion (HE) compounds, it is very common 
to establish the Kondo energy scale Tk from the linear co- 
efficient of specific heat 7 through Rajan's formula Tk= 
7tJR/3jo derived for the degenerate {2J + 1 > 2) Kondo 
modeli where J is the total angular momentum. In pre- 
vious studies of the specific heat of the HE compounds 
CesX (X=In, Sn)^''^ which crystallize in the CusAu cubic 
structure, this formula was used to determine the Kondo 
temperature, which was found to be 4.8 K for Cesin and 
16.7 K for CegSn. The crystal electric field (GEE) exci- 
tation energy was estimated to be Tcef=(!>^ K- 

Most HE compounds reside close to a quantum criti- 
cal point (QCP)where antiferromagnetic (AEM) or ferro- 
magnetic (EM) correlations are present. This makes the 
previous analysis inappropriate in so far as it assumes 
that the magnetic correlations do not contribute to 7. 
Indeed, the Wilson ratios (7r^i?xo/3Cj7o) which were de- 
termined previously for Gesin and GesSn are 11.5 and 7.0 
respectivelySJ^ indicating that ferromagnetic correlations 
dominate the susceptibility. 

Inelastic neutron scattering(INS) experiments on sin- 
gle crystals of compounds that are close to a QCP, such 
as CeRu2Si2'* or GeNi2Ge2- exhibit two classes of ex- 
citations. At most Q in the Brillouin zone, the scat- 
tering has the characteristic Kondo energy dependence 
and is Q-independent or only weakly Q-dependent. Sim- 
ilar behavior is observed in intermediate valence com- 
pounds for which it is clear that the behavior of the low 
temperature susceptibility, specific heat and INS spectra 
are close to the Kondo impurity prediction, as though 
the onset of lattice coherence has only a minor effect on 
these measurements^ii. Near the QGP, however, large 



Q-dependent scattering is observed with maximum in- 
tensity at the critical wavevector Qc {Qc = for EM 
and Qc = Qn for AEM) where ordering occurs in the 
nearby magnetic state. This scattering represents the 
short range order. It is dynamic and critically slows 
down, or softens, as the QCP is approached by lowering 
the temperature or changing a control parameter. These 
fiuctuations affect the specific heat and can result in non- 
Eermi liquid behavior. 

Hence INS in single crystals can separate the Kondo 
behavior from the contributions due to magnetic corre- 
lations. Since the spectral weight in the magnetic cor- 
relations is typically small, INS in polycrystals will be 
dominated by the Q-independent Kondo scattering. INS 
can also be used to directly determine the GEE excita- 
tions. Under these circumstance, INS provides a better 
way to determine Tk and Eqef than through analysis 
of the specific heat. In this paper, we employ INS to 
determine both Tk and Eqef- We have re- measured 
the magnetic susceptibility, and have extended the spe- 
cific heat measurement, which in the previous report was 
measured down to 1.8 K in zero applied magnetic field, 
to T = 400 mK and B = We have also measured the 
low temperature magnetization to 13 T. 

In the CeaX compounds, the Ce ions sit on the 
face centers of the cubic lattice and are subject to a 
crystalline electric field (GEE) of tetragonal symmetry. 
In this case, the Hamiltonian is described as: 



HcF = B^^Ol + BlOl 



where and 0™ are the crystal field parameters 

and Steven operators, respectively. The sixfold de- 
generate Af^, J=5/2 state splits into three doublets. 



2 



Diagonalizing the Hamiltonian. the atomic wave func- 
tions are given by^^ 

^rj\± 5/2 > +^1 -7721 T 3/2 > 

rf ) = v/l^l ± 5/2 > -n\ T 3/2 > 
re = I ± 1/2 > 

Depending on the admixture of the J=5/2 and 3/2 
states, the inelastic neutron scattering spectra will ex- 
hibit one or two inelastic excitations. Low energy trans- 
fer quasielastic scattering will also be observed if the in- 
strumental resolution is adequate. From the INS spectra, 
the crystal field energies and wavefunctions can be deter- 
mined from the ampitudes and energies of the excitations. 
The Kondo effect, which arises from the hybridization of 
the 4 /-electron with the conduction electrons, broad- 
ens the peak line-widths proportional to ksTx- The 
quasielastic scattering peak width Tqe can be equated 
to the Kondo energy UbTk of the ground state doublet. 

In what follows, we will use the CEF parameters and 
the Kondo energies derived from the neutron scattering 
to calculate the Kondo contribution to the specific heat, 
susceptibility and magnetization. All the Kondo calcu- 
lations utilizediiiSiiii employ the same Bethe-Ansatz cal- 
culation, making intercomparison possible. 

II. EXPERIMENT 

All samples were prepared by arc melting in an ultra- 
high-purity argon atmosphere. After melting the sam- 
ples were sealed under vacuum and annealed at SOO^C 
for 2 weeks and cooled slowly to room temperature. The 
magnetization was measured in a 14 T Quantum Design 
Vibrating Sample Magnetometer at the National High 
Magnetic Field Laboratory (NHMFL) at Los Alamos Na- 
tional Laboratory. The specific heat was measured in a 
Quantum Design PPMS system. The magnetic suscep- 
tibility measurements were performed in a commercial 
superconducting quantum interference device (SQUID) 
magnetometer . 

We performed inelastic neutron scattering on a 29 
gram sample of Cesin and a 37 gram sample of CeaSn us- 
ing the high energy transfer chopper spectrometer (HET) 
at ISIS (at the Rutherford Appleton Laboratory) and 
the low resolution medium energy chopper spectrome- 
ter (LRMECS) at IPNS (at Argonne National Labora- 
tory). For CcaSn, the quasi-elastic neutron spectrometer 
(QENS, at IPNS) was also used to measure the low en- 
ergy scattering. To increase the dynamic range of the 
INS spectrum, a variety of incident energies {Ei = 15 
meV, 35 meV, 60 meV, 100 meV for HET and 35 meV for 
LRMECS) and temperatures ( 4.5 K, 100 K, 150 K, 200 
K and 250 K for HET; lOK, lOOK, 150K for LRMECS) 
were employed. The HET data have been normalized to 
vanadium to establish the absolute value. All the data 
have been corrected for absorption ( which is very ob- 



vious for Ceain case), total scattering cross section, and 
sample mass. 

For the HET data, the low Q data were obtained from 
averaging the low angle detectors with angles ranging 
from 11.5 degrees to 26.5 degrees. The high Q data were 
obtained from the high angle detector bank at an angle 
136 degrees. For the LRMECS experiment, the low Q 
data were obtained by averaging over the low angle detec- 
tors with average angle equal to 13 degrees; and the high 
Q data were obtained from high angle detectors where the 
average angle was 87 degrees. The QENS data were col- 
lected at 7 K. This inverse geometry spectrometer has 19 
detector banks with Q from Q = 0.36 to Q = 2.52, each 
with a slightly different final energy {Ef from 2.82 meV 
to 3.36 meV). For every fixed Q, we removed the Ce 4/ 
form factor to obtain a spectrum representing the (5 = 
scattering and then summed all 19 spectra together to 
obtain a total S{Q = 0, /S.E) spectrum. In order to com- 
pare the QENS spectrum with the spectra from the direct 
geometry spectrometers HET and LRMECS, we multi- 
plied the QENS spectrum S{Q = 0, /\E) by the 4/ form 
factor apppropriate for HET at Ei = 35 meV. 

To subtract the nonmagnetic (background, single 
phonon, and multiple phonon) contributions, we mea- 
sured the non-magnetic counterpart compounds Laain 
and LaaSn. For the specific heat, we obtained the mag- 
netic contribution by direct subtraction, i.e. Cmag — 
C(Ce) - C{La). For the INS data we used Laain 
and LaaSn to determine the scaling of the nonmag- 
netic scattering between low Q and high Q as h{AE) — 
S{La, LQ)/ S{La, HQ). Using this factor we scaled the 
high Q data (where the nonmagnetic scattering domi- 
nates) to the low Q data (where the magnetic scattering 
dominates) in Ce compounds to determine the nonmag- 
netic contributioni^rii. 



III. RESULTS AND DISCUSSION 

Fig. 1(a) and (c) directly compare the low Q INS spec- 
tra of CeaIn and Laain; the data were collected on HET 
with incident energy Ei = 15 meV (a) and 60 meV (c) at 
4.5 K. Fig. 1(b) compares the low Q and high Q data for 
Ceain collected on LRMECS with an incident energy E.i= 
35 meV at 10 K. The low Q data for CeaSn and LaaSn, 
which were collected on HET, are compared in Fig. (d), 
(e) and (f) where the incident energies are Ei=15 meV 
(d). 35 meV (e) and 60 meV (f) at 4.5 K. In these spec- 
tra, two excited energy levels, corresponding to crystal 
field excitations, are observed for both CeaIn and CcaSn. 
The spectra (a) and (d), which compare the Ce3ln(Sn) 
and La3ln(Sn) scattering at low energy transfer {AE < 
9 meV), exhibit obvious quasielastic scattering which as 
mentioned above arises from Kondo scattering. 

The magnetic contribution Smag to the scattering 
of CeaIn, obtained using the method described above, 
is shown in Fig. 2. The solid lines represent a fit 
to the CEF model. Since the inelastic peaks are 
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FIG. 1: Inelastic neutron scattering spectra for Cesin and 
CeaSn together with that of their nonmagnetic counterpart 
compounds obtained from HET and LRMECS. The data col- 
lected on HET are at 4.5 K and on LRMECS are at 10 K. 
(a) Ei=15 meV, (b) Ei=35 meV, and (c) Ei=60 meV spectra 
of Cesin and Laaln. (d) Ei=15 meV, (e) Ei=35 meV and (f) 
Ei=60 meV spectra for CesSn and LasSn. All spectra are for 
low Q except in (b) where a high-Q spectrum is included for 
comparison. 

relatively broad, the line widths Ti are taken to be fi- 
nite. In this case, the magnetic scattering is described as: 

X"(Q, AE) ^ ^x,{T)AE{^)/[{AE - E^f + Tf] 

Here i=0,l,2, Eq = Q corresponds to the quasielas- 
tic scattering, and P{Q) is the Ce 4/ form factor. 
The CEF model fitting was performed simultaneously 
on six data sets at three different incident energies 
(£'i=15 meV, 35 meV and 60 meV) and at two different 
temperatures (4.5 K and 150 K). Fig. 2(a)-(d) are the 
data collected on HET. In Fig. 2(f) the LRMECS data 
are displayed for comparison. The resulting CEF fitting 
parameters are shown in Table I. The ground state is the 
Fy^'' doublet, the first excited state is the T^^^ doublelji^ 
at the energy 13.2 meV, and the second excited state is 
the Fg doublet at the energy 44.8 meV. The quasielastic 
line width Tqe = Fq = 1.49 meV, implies that Tr- = 
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TQE/kB = 17 K. 



TABLE I: CEF model fitting parameters for Ceain and CeaSn. 





Ceain 


CeaSn 


B2°(meV) 


-2.203±0.015 


-1.660±0.017 


BlijaeV) 


0.066±0.001 


0.038±0.0009 


B|(meV) 


-0.154±0.004 


-0.263±0.003 


V 






El (meV) 


13.2 


18.5 


E2 (meV) 


44.8 


36.1 


Tqe (meV) 


1.49±0.07 


3.52±0.16 


Fi (meV) 


5.98±0.07 


9.37±0.038 


Fa (meV) 


2.06±0.37 


6.28±0.38 


x' 


2.4057 


2.1528 


A (mole-Ce/emu) 


62 


85 



In Fig. 3(a)-(e) we display the magnetic contribution 
to the CesSn scattering collected from HET at 4.5 and 
100 K and at three incident energies (15, 35, and 60 meV). 
Data from QENS at 7 K (Fig. 3(f)) are included for 
comparison. The CEF fits are also included (solid lines) ; 
as for the Ce3ln case, the fits were performed simultane- 
ously on six different spectra at different incident energies 
and temperatures. The intensity and form factor of the 
QENS data have been adjusted to that of the HET spec- 
tra at £"^=35 meV (spectra (b)) to make a direct com- 
parison. The fitting parameters yield a similar crystal 
field scheme as for Ce^ln: the Tj^^ doublet is the ground 

(2) 

state, Fy ' is the first excited state with energy 18.5 nieV, 
and the second excited state is the Fg doublet at the en- 
ergy 36.1 meV. The Kondo temperature = 40 K is 
higher than for Ccsln, and the excited state linewidths 
are broader, reflecting stronger 4/-conduction hybridiza- 
tion. 

Due to the large CEF excitation energies, the low tem- 
perature behavior of the magnetic specific heat should be 
dominated by the Tj^^ doublet ground state. This is con- 
firmed by the fact that the magnetic entropy (Fig. 4(a) 
inset) reaches Rln2 near 20 K but only reaches Rln4 near 
70 K. For a doublet ground state, the Kondo model pre- 
dicts 7o — -kR/QTx for the linear coefficient of specific 
heat^. In previous results for Cesin a value Tk = 4.8 K 
was deduced using this formula^. In addition, the specific 
heat coefficient C/T showed a peak near 2 K whose exis- 
tence was somewhat uncertain since the lowest measured 
temperature was only 1.8 K. We have extended the spe- 
cific heat measurement down to 400 mK. In Fig. 4(a) we 
plot Cmag/T and find a peak at T = 2.6 K. Comparison 
of the data to the prediction 7^ [T) of the Kondo model 
which is calculated using the value Tk = ^QE/ks = ^IK 
deduced from our neutron data, shows that the Kondo 
prediction is much smaller than the experimental value; 
indeed, is only half of 7o^^ (Table II). Given the 
large Wilson ratio reported earlier^, the obvious expla- 
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FIG. 2: Magnetic contribution Smag to the inelastic neutron 
scattering spectra of Cesin for data taken on HET at T=4.5 
K and 150 K with different incident energies Ei = 15 meV, 35 
meV, and 60 meV and taken on LRMECS at T=100 K with 
Ei=35 meV. The solid lines represent the quasielastic and 
crystal field contributions obtained from least squares fitting 
as described in the text. 

nation is that ferromagnetic (FM) fluctuations dominate 
the low temperature specific heat, increasing the specific 
heat above the Kondo value and giving rise to the peak 
at 2.6 K representing the onset of short range FM order. 

We next consider the high temperature susceptibility, 
comparing the measured value to the value calculated 
from the crystal field parameters of Table I in the inset 
of Fig. 4(b). A molecular-field A = 62 mole-Ce/emu 
has been added to compensate the reduction of the sus- 
ceptibility at high temperature due to the Kondo effect 
= ^/x'^^^ + A). At high temperatures, when 
the crystal field states are excited, the effective Kondo 
temperature T^-^ is larger than the Kondo temperature 
of the ground state doublet. The molecular field con- 
stant is related to the effective Kondo temperature via 
A = where C5/2 is the free ion Curie constant 

for cerium. This relation gives = 77 K, which value 
is essentially equal to the width Fi of the first excited 
level seen in the neutron scattering (Fig. 2 and table I). 

At low temperatures, there should be three contribu- 
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FIG. 3: (a)-(f) Magnetic contribution Smag to the inelastic 
neutron scattering spectra for CesSn. The temperatures and 
incident energies are given in the plot. The solid lines repre- 
sent the CEF model, (f): the magnetic contribution to the 
INS spectra collected from QENS at 7 K. The solid line in (f) 
is the CEF model fit for the Ei = 35 meV spectra. 



tions to x{T), as well as to M{H) and Cmag- one from 
the Kondo single ion impurity physics of the ground state 
doublet, one from the FM fiuctuations, and one from the 
excitation of higher lying crystal field states. To carry 
out such an analysis, we note first that in the CusAu 
crystal structure, the tetragonal crystal field axis (i.e. 
the z-axis for the doublet wave functions) points per- 
pendicular to the face containing any given face-centered 
cerium atom; hence there are three orthogonal tetrago- 
nal axes in the unit cell. When applying a magnetic field 
in a polycrystalline sample, the field will point along the 
tetragonal axis for i of the cerium atoms but orthog- 
onal to the tetragonal axis (in the x ~ y plane) for | 
of the atoms. The effective low temperature Curie con- 
stant is then Cg// = |C|j.^ -I- ^C^jf, where Cl^fj = 
A^((/g|y''/iB)^i(i -I- 1)/3A:b. This is the form for a pseudo 
spin i doublet where the CEF physics is absorbed into 

the effective gr-factor. Here ffe//' ~ T Jz(x) > where 
< Jz{x) > is the matrix element of the angular moment 
component along the z{x) axis. From the CEF mixing 
parameter ry, we determine C^ff to be 0.48 emu-K/mole- 
Ce for Cesin and 0.41 emu-K/mole-Ce for CesSn. (Table 
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FIG. 4: (a) The magnetic contribution to the specific heat 
Cmag/T versus T for Cesln. The sohd Une is the Kondo pre- 
diction 7^(r) calculated using Tk = Tqe/Ub = 17 K. The 
inset is the magnetic entropy of Cesln. (b) The magnetic 
susceptibility x{T) for Cealn. The solid line is the Kondo 
prediction {T) calculated using Tk = 17 K and the low 
temperature Curie constant determined as described in the 
text. The dashed line gives the sum of the Kondo and crystal 
field contributions ^^^^bf+k = + [^cef _ ^lt^^j ^j^^ 

inset is the inverse susceptibility together with the calculated 
susceptibility (solid line) 1/x^^ = l/x'^^^ + -^ obtained using 
the CEF fitting parameters in table I. (c) The magnetization 
for Cesln. The solid line (B) is the Kondo calculation cal- 
culation using Tk ~ 17 K. The dashed line is the contribution 
from the ferromagnetic fluctuations. 

II). 

To sort the low temperature susceptibility into Kondo, 
FM, and CEF contributions, we note that since the first 
CEF excited level is at 152 K, at sufficiently low tem- 
peratures tlie Zeeman splitting of the ry^"* ground state 
doublet will obey a Curie law Xcurie — ^eff -^^^ 
the Kondo effect, this Curie behavior will be replaced by 
the Kondo contribution , which we calculate using the 
same Curie constant Cf^-y and using Tk = 17 K (solid 
line. Fig. 4(b)). The susceptibility from the combination 
of the ground state Kondo and the excited crystal fields 
will then be of the form x^^^+^ = + {x^^'' -Xclr^e) 
where we subtract Xchirie avoid double counting the 
ground state contribution. As for the specific heat, the 
resulting T^C'-^F-i-if (hashed line in Fig. 4(b)) is much 
smaller than the experimental value at T < 20 K. The 



FIG. 5: (a) the magnetic specific heat Cmag/T versus T at 
B=0 T (open circle) and B=9 T (solid triangle) for CesSn. 
The solid line is the Kondo contribution 7^(r) for Tk = 40 
K. The inset is the magnetic entropy, (b) Magnetic suscepti- 
bility x{T) for CeaSn. The solid line is the Kondo contribu- 
tion x^(T) calculated with Tk = VqE/kB = 40 K and the 
dashed line is the sum of the Kondo and CEF contributions 

^CEF+K = + (^CBF _ ^LT^.^y ^^^^^ ^j^^ -^^^^^^^ 

susceptibility together with the value 1/x^^ = ^/x^^^ + 
(solid line) calculated from CEF fitting parameters in table I. 
(c) Cmag at _B=0 T (opeu circle) and B=9 T (solid triangle) 
for CesSn. The thin solid Une is C^(T) calculated with Tk = 
40 K. 

excess can be viewed as the contribution from the ferro- 
magnetic fluctuations. Taking the latter to be equal to 
the difference x^""^ - X^'^''^^{T), the FM contribution 
is seen to increase below 10 K in a manner characteristic 
of ferromagnetic short range order. 

In Fig. 4(c) we exhibit the magnetization as mea- 
sured up to 13 Tesla at T=2 K. Based on Hewson's cal- 
culation of Kondo modeliS, we can estimate the Kondo 
contribution to the magnetization. Since the effective g- 
factors differ in the z and x — y directions, we calculate 
= iM^(z) + |M^(xy). The result is plotted as 
a solid line in Fig. 4(c). After subtracting the Kondo 
contribution, we obtain the contribution from the FM 
correlations (dashed line). This saturates at a relatively 
smaU field B ~ 2.5 tesla with M^"* = 0.095 hb, which is 
10 percent of the saturation value 1.0 expected based 
on the effective g-factors. 

In the same way as for Ccsln, we calculate x^{T), 
X^^^(r), X^^^+^(r), 7^(r) and C^(r) for CeaSn, 
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TABLE II: Kondo single ion model calculation for Cesin and CeaSn. 



- 7 ^ , T ^ ]\,rsat ^LT ^ ,K { emu \ / J \ ^e^P ( emu \ J 
< '^z > < ■Jx > MqeF '^efS XO [T^;:^T^:rc^ ) \moleCeK^> XoAK\ mole- Ce ) 'OAK \ moleCeKi ■ 

Cegln 2.0344 0.7171 0.991 0.4765 0.0076 0.256 0.064 0.467 

CegSn 1.6684 0.9074 0.994 0.4086 0.0033 0.109 0.018 0.221 



comparing to the measured data in Fig. 5. The high 
temperature susceptibihty (Fig. 5(b) inset) can again 
be fit with the sum of the CEF contribution calculated 
using the parameters of Table I and a molecular field con- 
tribution (solid line). The value A = 85 mole-Ce/emu of 
molecular field constant implies an effective Kondo tem- 
perature at high temperature = 105 K, which again 
is essentially equal to the linewidth 9.4 meV of the first 
excited state seen in the neutron scattering (Table I). 

The solid lines in Fig. 5 represent the Kondo ground 
state doublet contributions. In Fig. 5(b), the dashed line 
is x'-^^^+^(r). The excess due to the FM correlations 
has a much smaller magnitude (^ 0.005 emu/mol-Ce) 
than for Cesin where the FM contribution is of order 
0.05 emu/mol-Ce. A similar statement holds for the FM 
contribution to the specific heat coefficient, which is of or- 
der 0.4 J/mol-Ce-K2 for Cegln but only 0.1 J/mol-Ce-K^ 
for CesSn (Figs. 4(a) and 5(b)). Hence the FM enhance- 
ment is smaller in CcaSn than in Cesin, consistent with 
the larger value of Tk- 

In order to better understand these compounds, we 
measured the specific heat of Cesin under different ap- 
plied fields (B = T, 1 T, 3 T, 6 T and 9 T). The resuhs 
for the magnetic contribution C'mag are shown in Fig. 6. 
The low temperature peak in Cmag moves to higher tem- 
perature when the field is increased. Since the peak in 
the Kondo contribution to Cmag is expected to increase 
with field, we plot in Fig. 6(a)-(e) the Kondo contribu- 
tion (B) calculated for different applied fields using 
the theoretical results of Sacramento and Schlottmann^^. 
In calculating (B), we again account for the different 
effective (/-factors in the z and x — y directions. The 
results indicate that the Kondo contribution is not ex- 
pected to alter significantly in applied fields of order 9 
T, essentially because geff^^sB < kgTK for these fields. 
This makes it clear that the peak does not arise from the 
Kondo scattering but must be due to the FM fluctua- 
tions. 

To quantify the FM contribution, we again assume 
that the measured magnetic specific heat is the sum of 
the ground state doublet Kondo contribution {B), the 
FM contribution C^^' , and a contribution C^^^ due to 
the excitation of higher lying CEF states. Since the FM 
fluctuations appear to only contribute to the susceptibil- 
ity below 10 K (Fig. 4 b) we assume that the excess Cmag 
- C^{B) observed for T > 10 K is primarily due to CEF 
excitations. Given the large linewidths of the CEF exci- 
tations seen in the neutron scattering, and concomitant 
large effective T^^ at high temperature, this contribu- 
tion to the specific heat is much broader as a function 



of temperature than would be the case for a simple CEF 
Schottky anomaly. For simplicity, we approximate this 
CEF contribution as linear in temperature, with slope 
equal to that observed in the range 8-15 K, and we as- 
sume that since the CEF excitation energy is large, this 
contribution will be unaffected by fields of order 9 T. We 
approximate the FM contribution C^^^ as a Gaussian, 
centered at a temperature that increases with field. The 
three contributions, Kondo, CEF, and FM, are plotted 
at the different fields in Figs. 6(a)-(e). The solid lines, 
which represent the sum of all three contributions, fit the 
data very well at all fields. 

We plot the Gaussian peak temperature in Fig. 6(f), 
where it is seen to grow linearly with field. This suggests 
Zeeman splitting, where at zero field the splitting arises 
from the internal field in the regions of FM short range 
order, and where the applied field increases the split- 
ting. To determine the internal field Bint, we calculate 
the Schottky anomaly C^'^^"*'^'^^ {Bmt) expected due to 
Zeeman splitting of a doublet with the same effective g- 
factors as we have obtained from the neutron fits; we then 
adjust Bint until the peak temperature of the Schottky 
anomaly is the same as that of the Gaussian peak tem- 
perature for 5 = 0. This gives Bint = 9.5 T. We then 
calculate C^'''^°^*'^'^{Bint + Bapp) to determine peak po- 
sition of the Schottky anomaly in an applied field Bapp. 
As can be seen in Fig. 6 (f), the Gaussian peak temper- 
atures track the expected Zeeman splitting very closely. 
On the other hand, the temperature dependence of the 
Schottky specific heat calculated in this manner is con- 
siderably broader than the Gaussian contributions 
that are plotted in Fig. 5. This means that, while the 
contribution of the FM short range order to the specific 
heat is not of Schottky form, the increase of the Gaus- 
sian peak position is the same as the Zeeman splitting 
expected for a total field Bint + Bapp given the effective 
g-factors. 

The entropy of the Gaussian contribution is about 15 
% of Rln2 for all fields. This corresponds to the estimate 
obtained from the magnetization M (B) where the satu- 
ration value of the FM contribution is about 10% of the 
value 1.0 /is expected for the Tj^'^ ground state doublet 
for the measured value of rj. Hence, the enhancement of 
Xo and 7o arises from magnetic fluctuations which involve 
10-15 % of the 4/ electron degrees of freedom. 

In Fig. 5(c), we compare the magnetic specific heat 
Cmag at zcro field and i3 = 9 T for CeaSn. The solid 
line is the Kondo contribution C^ . The specific heat 
does not change with field for i? < 9 T. The most likely 
explanation of this is that, as discussed above, the FM 
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FIG. 6: (a), (b), (c), (d) and (e): Cmag of Ceain in different 
applied magnetic fields. The solid line sums the three con- 
tributions (Kondo, CEF, and FM fluctuations) shown in the 
plot, (f): The peak position of Cmag, of the Gaussian con- 
tribution C^^^ due to the FM fluctuations, and the expected 
peak position in the Schottky anomaly 0^"''°^^''^ {Bi„t + Bapp) 
due to Zeeman splitting in the presence of an internal fleld 

Bint ■ 

correlations make a smaller contribution than in Ccsln. 
The excess specific heat Cmag — seen for T > 6 K is 
presumably due to the CEF contribution, which should 
be even broader in temperature in CcsSn than in Ceain 
due to the larger Kondo temperature. 

We have demonstrated that the large 7 observed in 
Cesin arises more from ferromagnetic correlations than 
from the single ion Kondo physics. This reflects the fact 
that the system is close to a ferromagnetic quantum criti- 
cal point. In a Q-resolved INS experiment, the ferromag- 
netic correlations should show up in the vicinity of Q —0 
riding on a background of Q independent Kondo scatter- 
ing. We can estimate that these FM correlations will have 
10-15% of the total spectral weight in Q-spa.ce. Since the 
large 7 and the proximity to the QCP occurs when the 
Kondo temperature Tk = 17 K is fairly large, it is also 



reasonable to believe that when an appropriate control 
parameter (e.g. alloying parameter x in Ce3_j;Laa;In) 
drives this system to the QCP, the Kondo temperature 
Tk will remain finite, as expected for example for a spin 
density wave type QCP. 

We have observed strong FM fluctuations in the related 
compound Prsin, which is an antiferromagnet below 12 
Ki^. A possibility for this behavior is that AFM inter- 
actions between rare earth atoms on the face centers of 
the CusAu structure are frustrated. If, for example, the 
atoms at (1/2 1/2 0) and (1/2 1/2) are aligned anti- 
ferromagnetically, the atom at (0 1/2 1/2) will be free 
to point to any direction. Ferromagnetic next-nearest- 
neighbor interactions could then stabilize ferromagnetism 
on this sublattice^^. In any case, the FM correlations ap- 
pear to be generic to this crystal structure. 

In conclusion, we have used inelastic neutron scattering 
to determine the crystalline electric field (CEF) splitting 
and Kondo energy scale in Cesin and CeaSn. For both 
compounds the crystal field excitation energy is large. 
For Ccsln we have separated the magnetization M{B), 
susceptibility x{T^) ^-^id specific heat Cmag into contribu- 
tions from the Kondo effect, from the CEF, and from FM 
fluctuations. The simplified model calculation for Ceain 
shows that the FM correlations make a 15 % contribution 
to the doublet ground state entropy and that the large 
7 arises mostly from the FM correlations. This suggests 
Ceain is close to a quantum critical point (QCP). The 
Kondo temperature Tk is expected to remain finite at 
the QCP, as occurs for a spin density wave type QCP. 
INS experiments in single crystals of these compounds 
would be very interesting. 
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